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Abstract 

In this paper, we study two different nonhnear interpolating paths in adiabatic evo- 
• lution algorithms for solving a particular class of quantum search problems where both 

the initial and final Hamiltonian are one-dimensional projector Hamiltonians on the cor- 
responding ground state. If the overlap between the initial state and final state of the 
. quantum system is not equal to zero, both of these models can provide a constant time 

speedup over the usual adiabatic algorithms by increasing some another corresponding 



"complexity" . But when the initial state has a zero overlap with the solution state in the 
problem, the second model leads to an infinite time complexity of the algorithm for what- 
ever interpolating functions being applied while the first one can still provide a constant 
running time. However, inspired by a related reference, a variant of the first model can be 
constructed which also fails for the problem when the overlap is exactly equal to zero if 
we want to make up the "intrinsic" fault of the second model — an increase in energy. Two 
concrete theorems are given to serve as explanations why neither of these two models can 
improve the usual adiabatic evolution algorithms for the phenomenon above. These just 
tell us what should be noted when using certain nonlinear evolution paths in adiabatic 
quantum algorithms for some special kind of problems. 
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1. Introduction 



Quantum adiabatic evolution [1] has been developed as an alternative to the traditional 
circuit model of quantum computation and the equivalence between these two models has been 
proven recently |2], |3]. But adiabatic quantum computation has demonstrated some of its 
own advantages, especially like the robustness against different kinds of perturbations, such as 
decoherence [3], and unitary control errors [5J. 

The main idea behind an adiabatic quantum algorithm can be summarized as follows: 
prepare a quantum system in the ground state of a simple initial Hamiltonian, then slowly 
change this Hamiltonian to a final Hamiltonian whose ground state encodes the solution to 
the problem that will be solved. More precisely, assume the Hamiltonian H{t) of a quantum 
system can fully characterize the time evolution of the system. Let ly^it)) be the state of the 
system which is evolving according to the Schrodinger equation 

^J^\^{t)) = H{t)\^{t)). (1) 

Define \E{0,t)) and |-E'(l,t)) as the ground state and first excited state of the system Hamil- 
tonian H(t) with corresponding eigenvalues E{0,t) and E{l,t), respectively. The quantum 
adiabatic theorem [B] tells that if we prepare the system in the ground state of the initial 
Hamiltonian Hi and let it evolve slowly along some path to the final Hamiltonian Hf, then 

|(E(0,r)|^(T))r>l-e2(0<5«l), (2) 

provided that 



mm 



SI. 

in which 



max 

0<t<T 



{E{l,t)\^\E{0,t)) 



(4) 



is a measurement of the evolving rate of the Hamiltonian, and 



S^in=min^[E{l,t)-E{0,t)] (5) 

is the minimum gap between the two lowest eigenvalues. In general, the required run time 
T for a typical adiabatic algorithm will be mainly determined by 6min so long as Amax is 
polynomially bounded. For the sake of convenience, the time-dependent Hamiltonian H{t) can 
be reparametrized as 

H{s) = H{t/T), /or < s < 1. (6) 



Usually, the instantaneous Hamiltonian H{s) which connects the initial Hamiltonian Hi and 
final Hamiltonian Hf of the quantum system is a linear adiabatic evolution path 

His) = (1 - s)H, + sHf. (7) 

But this simple form for the interpolating Hamiltonian may sometimes have some limitations 
for some special problems, as we will show. 

In fact, there is no reason not to try paths involving terms that are not linear combinations 
of the initial and final Hamiltonians. The adiabatic algorithm will work taking any path H{s), 
as long as the adiabaticity condition is satisfied: the ground state at s = is the initial state 
and the ground state at s = 1 is the solution state. In Ref. the authors had proposed the 
following general form of time-dependent Hamiltonian, 

H{s) = (1 - s)Hi + sHf + s{l - s)He, (8) 

in which H^ is an extra piece of the Hamiltonian that is turned off at the beginning and end 
of the adiabatic evolution. It was found this nonlinear evolution can change the performance 
of the algorithm from unsuccessful to successful when solving some special problems, such as 
those discussed in Ref. [8]. As a result, when applying this method to the quantum search 
problem [H], by a suitable choice of the time- dependent Hamiltonian, it is possible to do the 
calculation in constant time, independent of the number of items in the database. However, 
in this case, the original time-complexity is replaced by the complexity of implementing the 
driving Hamiltonian He [TO] . 

In Ref. [TT], a general class of models for the time- dependent Hamiltonian was proposed, 
which has a form like: 

H{s) = f{s)H, + g{s)Hj, (9) 
where f{s) and g{s) are arbitrary functions of time, subject to the boundary conditions 

/(O) = ^7(1) = 1, f{l) = g{0) = 0. (10) 

When performing unstructured adiabatic quantum search of a database with items, it was 
found that by increasing the lowest eigenvalue of the time- dependent Hamiltonian temporarily 
to a maximum of oc \/N, it was possible to do the calculation in constant time as well. 

In this paper, we study a special class of quantum search problems where both the initial 
and final Hamiltonian are one-dimensional projector Hamiltonians on the corresponding ground 



state by using the above two different approaches. When the initial state has a nonzero overlap 
with the final state, both of these approaches can provide a constant running time by increasing 
some another corresponding "complexity" . However, if the overlap is exactly equal to zero, we 
find that only the adding a driving Hamiltonian approach whose form is given by Eq. can 
still work while the other approach fails no matter what the interpolating functions chosen for 
f{s) and g{s). But at the same time, if we want to eliminate the pitfall of the general class 
of models for the time-dependent Hamiltonian — an increase in energy for the system, we find 
that a variant of the driving Hamiltonian model can be built which also lead to an infinite time 
complexity for the problem when the initial state of the system is orthogonal to the solution 
state. These two facts just tell what we should note when using some nonlinear interpolating 
paths in adiabatic evolution algorithms for some special kind of problems. 

The organization of the current paper is as follows: in part 2, we define the problem that will 
be studied throughout the paper. The two different approaches for the algorithmic performance 
speedup of this problem as compared to a local adiabatic [i2l or partial adiabatic evolution 
algorithm [13] are discussed detailedly in part 3. Two concrete theorems are given to show 
why both of the variant version of the first model and the second model fails for the problem 
when the initial state is orthogonal to the final state in part 4. We summarize and draw some 
conclusions in the last part. 



For the purpose of illustrating our main idea, we discuss the adiabatic evolution for the 
problem studied in the current context which has a form for the initial and final Hamiltonian 

as 



It is easy to know that in fact our this defined problem can be seen as a kind of abstraction 
of quantum search problem. Now the problem is how long it will take us to evolve the initial 
state la) to the final state \f3) by applying a nonhnear global adiabatic evolution algorithm. In 
fact, if a simple linear adiabatic algorithm [13] is applied, it is easy to know the running time 
can be estimated as 



2. The problem 



Hi = I — \a) {a 



Hf = I-\[3){P\. 





(12) 



where 

a = {a\(3). (13) 

While a local adiabatic algorithm [12] or a partial adiabatic algorithm [13] is used for the same 
problem, it is not difficult to see that the time complexity is 

T2 = 0{\a\-'), (14) 

which provides a quadratic speedup over the global adiabatic evolution. But can we further 
reduce the time consumption for the problem? Inspired by the two different approaches intro- 
duced in the ffist section, we give the answer to this question in the positive side, as we will 
show in detail below. 

3. Two different models for the algorithmic speedup of 
the problem 

In this section, we show how the two different strategies introduced above can provide a 
speedup over the usual adiabatic algorithm such as global adiabatic, local adiabatic or partial 
adiabatic algorithms in subsection 1 and 2, respectively. Before we start, we introduce the 
following notations of a reduced two dimensional- Hilbert space to simplify the calculation of 
the two lowest eigenvalues of the Hamiltonians that will be shown: 

|l) = |a), |2) = i(|/3)-a|a)), (15) 

in which 

h = v^l - |a|2. (16) 

Thereby, we can get 

|«) = |1), |/3) = a|l) + 6|2). (17) 

3.1. Adding a driving Hamiltonian 

Consider a quantum Hamiltonian which has the form of Eq. ([8]) for solving the problem, in 
which Hf. is specified as: 

H,= \a){(3\ + \l3){a\. (18) 



Now the Hamiltonian H{s) in Eq. ([8]) has a matrix form 

/ sb'^ + 2s(l- s)Re(a) sb(l - s - a) . 
H{s) ={ ' ' ^ ^ I , (19) 

y sb{l - s-a*) 1 - sb^ 

by the notations being introduced at the beginning of this section, and the eigenvalues of it are 
easy to obtain: 

^ , ^ A- B ^ , ^ A + B 

Eo{s) = , E,{s) = , (20) 

in which we have used the following marks for simphcity: 

A = l + 2s{l- s)Re{a), (21) 



5 = y 1 - 4s(l - s)|i?e(a) + b'^ + s{l - s)[Im{a)]^ - s{l - s)}. (22) 
We therefore can get the energy gap of the system Hamiltonian: 



Ei{s)-Eo{s) = yi-4s(l-s)|i?e(a) + 62 + s(i_s)[/m(a)]2-s(l-s)|. (23) 
And the minimum energy gap reads 



5™„ = min [Ei{s) - Eo{s)] = J \a\^ - Re{a) - hlm{a)]^ + - (24) 

0<s<l V 4 4 

Thereby, when the overlap a is tending to zero or exactly equal to it, the time complexity 
for the algorithm tends to a constant. Simultaneously, it is obvious that the usual adiabatic 
algorithms such as linear adiabatic evolution etc. have a bad performance at this time. In 
fact, if a = holds, this will lead to the total failure of all the usual adiabatic algorithms — i.e., 
infinite time complexity, which can be seen as follows: 

sb'^ —sab 

H{s) = I , (25) 



-sa*b 1 — sb 



2 



SO the energy gap of it now reads: 



SLn = minVl - 4s(l - s)il - |a|2) = 0. (26) 

0<s<l 

But in the meanwhile, an obvious fact should be easy to observe: the added driving Hamiltonian 
in this model also takes additional "complexity" to implement compared with the widely studied 
adiabatic evolution procedure in the literature. However, this just implies more feasibilities can 
be achieved for the adiabatic evolution based quantum algorithms than its equivalent quantum 
algorithms based on the circuit model. 



3.2. Using a more general model of interplaiting path 

Now we turn to another approach to solve the problem. Consider the general interpo- 
lating form for the Hamiltonian H{s) in Eq. dH]), which has a matrix form in the reduced 
two-dimensional Hilbert space defined before: 

/ qb'^ —qab 

H{s) =\ ' ^ I , (27) 

y -ga*h f + g\aY 

and its two corresponding eigenvalues are therefore given by 



2 



E,{s) = ^^^^ )L^±_^ ii^^ (28) 



if + g) + ViT^W+W\ 



2 



So, we have 



when we set 



E^{s) = ^ ^ v^±_^ (29) 



X 

5min = min [Ei{s) - Eo{s)] = \a\ + -\a\, (30) 

0<s<l z 



f{s) = l-s + xs{l-s), g{s) = s + xs{l- s). (31) 



It is easy to see that if we choose a; = j^, the running time of the algorithm tends to a constant 
when the overlap a tends to zero. But simultaneously, it can be found that the ground state 
energy Eq{s) grows with a maximum value proportional to at the point s = | as well before 
it returns to zero at the end of the adiabatic evolution. In contrast, when the free parameter 
"x" disappears, this model has degenerated to the linear adiabatic evolution version, and the 
performance of it is very bad at this time (the overlap a tends to zero) as can be seen from 
the preceding subsection, although the ground state energy of the quantum system tends to a 
constant at the minimum time point. In fact, when the parameters "f(s)" &"g(s)" are specified 
as in Eq. ( I3T|) . we have the following even more accurate relation between this general model 
of interplaiting path and the usual linear adiabatic evolution: if the minimum gap between the 
ground state and the first excited state of the former one increases to an amount approximate to 
"x" larger than the second one, the corresponding ground state energy is enlarged to the same 
amount accordingly. And this is consistent with our intuition: adiabatic evolution performs 
slowly around the transition point where the energy gap is small, and if large energy is injected 
into the quantum system, it will be helpful for improving the performance of the adiabatic 
algorithm. 



What about the time complexity of this general class of model of adiabatic evolution algo- 
rithm for the current problem when the overlap a is exactly equal to zero? Obviously, Eq. ( 13T|) 
will still give an infinite time complexity as the usual adiabatic algorithms do no matter what 
value is taken for the free parameter "x". But what about the result if other functions for f{s) 
and g{s) are chosen? We will discuss this in the next section. 



4. The theorems 

In fact, we can prove the second approach always fails for the problem when the initial state 
has a zero overlap with the final state. 

Theorem 1 The second model of the interpolating Hamiltonian discussed in section 3.2 
will always cause an failure the adiabatic algorithm for whatever functions of f{s) and g{s) are 
selected when it is used to solve the problem above if the initial state is orthogonal to the final 
solution state. 

Proof By Eq. (125]) and Eq. (1291) . it is known that once the function /(s) and g{s) are 
chosen, the minimum energy gap between the lowest and second lowest eigenvalues is 

5^in = min - Eo{s)] = min |/(s) - ^(s)|. (32) 

0<s<l 0<s<l 

From the boundary conditions given in Eq. f lTOj) . a new function can be constructed as 

Fis) = fis)-g{s), (33) 

and it follows that 

F(0) = /(O) - g{0) = 1, F(l) = /(I) - g{l) = -1. (34) 

As functions f{s) and g{s) are both continuous and by "Intermediate Value Theorem" in 
mathematics, we immediately know that there must exist a point for which 

F{s) = 0, (35) 

holds for some < s < 1. This just implies an infinite time complexity for the adiabatic 
algorithm and the proof is completed . □ 
As we have already seen from section 3.2, the speed up of the adiabatic algorithm for the 
problem requires an increase in energy, at least temporarily. If we want to keep the ground 



state energy zero as stated in Ref. [TJ, we could use the following adiabatic Hamiltonian as 
shown in the next theorem, which can be seen as a variant of the first model discussed in section 
3.1. However, it also has the limitation that it fails for the problem when the overlap between 
the initial state and final state is exactly equal to zero. 

Theorem 2 The following adiabatic Hamiltonian form still fails for the problem defined 
above when the initial state has a zero overlap with the final state: 

H'is) = f{s)Hi + g{s)Hf - h{s)I, (36) 

in which f{s) and g{s) are arbitrary continuous functions, and h{s) is given by 

Ms) = /W + gW-l/M-9WI (3,) 

Proof Firstly, it can be easily verified that h{s) satisfies the correct boundary condition 
when adding an extra Hamiltonian model, the driving item only appears during the adiabatic 
evolution process: 

h{0) = h{l) = 0. (38) 
By the same method as before, we can write out the matrix form for the Hamiltonian H'[s) as 

— / g-h{s) \ 

H'{s) =1^'^ . (39) 

V f-h{s) J 

It follows that the minimum energy gap can also be easily obtained yielding 

6min = min |/(s) - g{s)\. (40) 

0<s<l 

By the same proof progress of theorem 1, we obtain the conclusion that the time complexity 
for this adiabatic algorithm is also infinite. □ 



5. Conclusion 

We have studied two different models of nonlinear paths in quantum adiabatic evolution 
for solving a special kind of quantum search problem in which both the initial and final Hamil- 
tonian are a one- dimensional projector Hamiltonian on the corresponding ground state. When 



the initial state has a nonzero overlap with the solution state, both of these approaches can 
provide a constant running time although the mechanisms behind them are different: in the 
first approach, the time complexity has been replaced by the "complexity" of implementing 
the driving Hamiltonian that is added while the ground state energy of the quantum system 
grows maximally proportional to a corresponding quantity temporarily before returning to zero 
at the end of the adiabatic evolution in the second model. But when the overlap between the 
initial state and the final state is exactly equal to zero, only the first approach can still work 
for the problem. A concrete theorem is proven to explain why the second approach fails for 
whatever interpolating functions are selected in that model. However, we have also shown that 
a variant version of the adding a driving Hamiltonian model which was mentioned in some 
related literature, will still lead to an infinite time complexity for the problem in the occasion 
that the overlap is equal to zero when we want to avoid an increase in energy temporarily. 
Using some nonlinear evolution paths in adiabatic algorithms may not always be helpful for 
certain special kind of problems, and we just should pay attention to this phenomenon when 
designing the adiabatic quantum algorithms. 
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